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Abstract 



> 

We exhibit a general procedure to purify any given ensemble by identifying an appropriate inter- 
action between the physical system S of the ensemble and the reference system K. We show that 
the interaction can be chosen in such a way to lead to a spatial separation of the pair S-K. As a 
consequence, one can use it to prepare at a distance different equivalent ensembles. The argument 
associates a physically precise procedure to the purely formal and fictitious process usually consid- 
■ ered in the literature. We conclude with an illuminating example taken from quantum computational 

theory. 

9*'- 

1 Introduction 

A statistical ensemble £ of physical systems 5 is characterized by a (finite, countable or continuous) set of 
positive numbers pi summing up to 1 and by a corresponding set of normalized vectors \tpi) of the Hilbert 
space H s associated to the system S, so that we will write £(pi, \tpi}) to represent it. The statistical 
operator ps (a trace-class, trace one, semipositive definite operator) associated to £(pi, |t/>t)) is defined 

H ' 

C3 1 as: 



P£=5>ihWil- (1.1) 



A point of great conceptual relevance which marks a radical difference between the classical and 
quantum cases is that, while in classical mechanics the assignment of the statistical operator p(r, p) 
uniquely identifies the ensemble, within quantum mechanics, as it is well known, the correspondence 
between statistical ensembles and statistical operators is infinitely many to one. 

With reference to this point, let us consider the set of all statistical ensembles of systems like the one 
under consideration. Such a set can be naturally endowed with an equivalence relation. 

Definition: We will say that two statistical ensembles £ and £* are equivalent, and we will 
write £ = £* , iff p £ — p £ *. 



* e-mail: bassi@ictp.trieste.it 
'e-mail: ghirardi@ts.infn.it 



1 



It is obvious that the just defined relation is reflexive, symmetric and transitive and that it leads to 

a decomposition of the set of all ensembles into disjoint equivalence classes. We will denote as [£} the 

equivalence class containing the ensemble £. 

Purification of an ensemble is a procedure by which one associates to the ensemble a pure state 
of an appropriately enlarged Hilbert space H s+K = H s ® TL K , where if is a reference system whose 

Hilbert space Tt K we assume to be infinite-dimensional for reasons which will become clear in a moment. 

The fundamental request on is that, by measuring an appropriate observable of K and confining 

attention to the system S alone, one can prepare the desired ensemble £(pi, |t/'*))- 

The first proof that, given two equivalent ensembles £{di, \ <pi)) and £{pj,\Xj))i one can fi n d two 

orthonormal sets {|^4i)} and {|-Bj)} oiH. K such that 

i*> = E \/*i&> ® = E vpj\xj) ® i^) (1.2) 

» i 

has been exhibited by Gisin Q . This result is particularly relevant since it is related to the request that 
no faster-than-light signals can be send between distant observers. 

Subsequently, Hughston et al. || have generalized the above result, providing a complete classification 
of equivalent ensembles: using the purification procedure, they have derived necessary and sufficient 
conditions for two ensembles to be equivalent. 

In the literature (see, e.g., ensemble purification is usually considered as a purely mathematical 
tool: one does not identify any dynamical mechanism which could be used to actually implement it, and 
the system K is considered a fictitious system without a direct physical significance. The aim of this 
paper is to exhibit a precise physical procedure in order to purify any ensemble by making the system S 
interact with a system K, in such a way that the desired pure state \^>) be actually produced. Then one 
can use it to prepare any desired ensemble of the equivalence class. 

2 Statistical ensembles and the purification process: the con- 
structive procedure 

As remarked above, it is our purpose to present a formal constructive mechanism to purify any given 
ensemble, showing at the same time how, by resorting to this procedure, one can use the obtained pure 
state to generate all ensembles of systems S equivalent to the one one has purified. The procedure is 
based on a formalism which parallels strictly the one proposed by von Neumann for implementing ideal 
measurement processes of the first kind, even though the system K, which plays a role analogous to the 
one of the measuring apparatus in his treatment, can very well be (and actually we will consider it to be) 
a microsystem. 

Our starting point is the consideration of an equivalence class [£] of ensembles of systems S. Within 
such a class there is the ensemble £ (dj, \(f>i)) which corresponds to the spectral decomposition of the 
associated statistical operator having the positive numbers di as eigenvalues and the \4>i) as the associated 
orthonormal eigenvectors. Such a decomposition is unique, apart from accidental degeneracies which, if 
they occur, can be disposed of as one wants, so that we will consider the eigenvectors \<f>i) as precisely 
assigned vectors. We assume that the index i runs from to n, without committing ourselves about the 
fact that n is finite or infinite and about the fact that the orthonormal set {|</>i)} be a complete set of 
7i s or not. 
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Let us consider now the orthonormal states \<pi) and let us assume that there exist a physical system 
K, whose associated Hilbert space Ti K is infinite dimensional, a state \ao) of Ti K and an interaction 
hamiltonian H s+K of TL S+K such that the S-K interaction lasting for a certain time interval T induces 
the following evolution: 

\(f>i) ® K) => \4>i) ® \a%), (a%\aj) =Sij, (2.1) 

where |oj) are statevectors belonging to 

In the next section we will exhibit a simple hamiltonian having such a property and leading also to an 
arbitrarily chosen separation in space of the systems S and K. We stress that we need TC K to be infinite 
dimensional if we want to be able to build a state \^t) which will allow us to prepare any ensemble 
whatsoever in the equivalence class under consideration by measurement procedures on systems fsT, since 
in any equivalence class there are always ensembles containing an infinite number of states. 

Given the ensemble £ (d,, \4>i)) we consider the state: 

n 
i=l 



and we let it evolve through the interval T. According to Eq. (2.1) and due to the linearity of the 
quantum evolution, we get 

n 

|*0>=H*T>=53V*l&>®|Oi>- (2 ' 3) 



We consider now an arbitrary complete orthonormal set j = 0, 1, oo, of TL and we complete 

(if necessary) the set {|ai)}, i = 0, 1, ...,n, to a set {|^4i)} by adding to it orthonormal states spanning 
the manifold oiTL K orthogonal to the one generated by the {|ai)} themselves. Obviously we have: 

oo 

l^)=E^I^)' i = 0,l,-.°o, (2.4) 
where Uij is a unitary matrix of 7i . From Eq.(5) we get: 

n oo 

z__0 j=0 

oo n 

= ElE^'l^ 8 ^ 

oo 

= £l&)®l&i>. (2-5) 



Note that = 1 implies: 



EfeK & ii^)te) = Eiife)ii 2 = i- (2-6) 

j,fc=0 3=0 



The states \xj) are not normalized, so that, putting \xj) = we have: 

oo 

l*T)=Ellfe)IIIXi)®^)- (2-7) 
j=o 

If we measure now an observable of the system K having a non-degenerate spectrum with \bj) as 
eigenvectors and we confine our attention to the resulting ensemble of systems S, we obtain the ensemble 
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£(lllx.j)l| 2 , |Xi))- No te that since both E™=o d il^)(^l and £j=o IIIXj)l| 2 |Xj>(Xjl are obtained by taking 
the partial trace on Tt K of I^tX^tI, they are equal and the corresponding ensembles belong to the same 



equivalence class. Thus we have proved that starting from the state (2.S) and choosing an observable 
having \bj) as eigenstates, we generate an ensemble which belongs to the same equivalence class of the 
original one. 

The relevant question we have to face now is the following: can all statistical ensembles belonging 
to the equivalence class of £(di, \4>i)) be obtained by properly choosing the observables of the system K 
we are going to measure? The answer is yes, as it is easily proved. To this purpose, let us consider 
an arbitrary ensemble £{pj, equivalent to £(di, \4>i)); we suppose that the index j runs from to 
N (> n), without excluding the case in which N is infinite. We know that the fact that the statistical 
operators associated to such ensembles are identical implies that the normalized states |tj) are linear 
combinations of the orthonormal states \<fii): 

n 

bj>=J>i|0i>, j = 0,l,...,N. (2.8) 

i=0 

We define now a rectangular matrix Vy having n + 1 rows and N + 1 columns by putting: 



V v = \lT b *> i = 0,l,...,n; j — 0, 1, N. (2.9) 



From the relation ^™=o di\<t>i) {fal = Yli^o Pj\ T j)( T j\> usm g Eq. (10) we immediately get: 



The above relation implies: 



JV 

*^2pjbjib* k = diS lk . (2-10) 

3=0 



N N i — - , — - 

= E \It. \lf k hiih * = 5ik - (2 - n) 



3=0 3=0 

We thus have n + 1 normalized and orthogonal vectors {w r }, r = 0,1, n of C N+1 , whose components 
are the row elements of the matrix V r j : 

w r = (V r0 ,V rl ,...V rN ), r = 0,l,...,n. (2.12) 

If N is finite, we pass from the vectors {w r } to new vectors {w r } of C°° by considering equal to zero 
the components of {wy} from N + 1 on. We then extend the set {w r } to a complete orthonormal set 
of C°°, by adding appropriately chosen normalized vectors {w s }, s = n + 1, .., oo. Correspondingly, the 
rectangular matrix V{j of Eq. ( |2.9| ) is transformed into an infinite square matrix, whose rows are the 
components of the vectors {w r }, for r = 0, 1, .., oo. Due to Eq. ( [2.1l| ) and the procedure we have followed, 
this infinite square matrix — which we keep calling — is unitary. 

Let us consider now an observable fi of TL K having a purely discrete and non degenerate spectrum 
with eigenvectors \Bj) = V^Ai); this implies that \Ai) = Vlj\Bj). Since Vij is unitary, we can repeat 
the previous procedure which amounts simply in replacing, in Eq. ( [2.3| ), the states |a^) = 1^) appearing 
there with their Fourier expansion in terms of the set {\Bj}}. Then Eq. ( |2.3| ) takes the form ( |2.7| ) where, 
according to the definition of \xj) given in Eq. ( f2~5| ) and of |r,) given in Eq. (^): 

n n 

\xj) = E v^Kj l^i) = y/pj E b a 1^') = 

i=0 i=0 

Ti\rjY (2.13) 



4 



This shows that |||xj)|| 2 = Pj and that normalizing \xj) we get the states \rj). Accordingly, Eq. (2.7) 
becomes: 

N 

\*t) =J2VPJ\Tj) ® (2.14) 

3=0 

so that measurement of Q, K reduces the state Y&t) to the desired ensemble £ (j>j, \tj)). Since £(pj, 
is an arbitrary ensemble belonging to the equivalence class [£ (di, \4>i))], this completes our proof. 

The now obtained result shows how, once one has prepared the pure state |^t), he has an immediate 
complete classification of all ensembles belonging to the equivalence class of £ (di, \(f)i)), an alternative 
way of deriving the nice result of || . 

Of course, within any equivalence class there are also mixtures which involve a continuous union of 
pure states i.e. 

£(p(A),|M) — PE = J d\p(X) \(/>\}{<f>\\, (2.15) 

with J dXp(X) — 1. To get such mixtures from the pure state we have, obviously, to measure with 
infinite precision an observable of TL K having a continuous spectrum. This is formally but not practically 
feasible. 

Concluding, if we can implement our "von Neumann-like ideal interaction scheme" we can perform 
the desired purification and then prepare any one of the ensembles in the equivalence class of £ (di, \ 4>i)) 
by performing an appropriate measurement on the system K. 

3 The appropriate hamiltonian for the desired purification 

To face our problem let us consider the following self-adjoint operator of 7i s+K : 

Hj = i \4>j){4>j\ ® [\a )(a 3 \ - \ aj ){a \], (3.1) 
and let us evaluate its powers. We have: 

//;" ' : = //, 

Hf n = l&X&l® [\ao)(a \ + \a 3 )( aj \}. (3.2) 
Let us consider now the operator exjp(-iujHjT): 

expi-iuHjT) = cosiuHjT) - i sm(wHjT). (3.3) 
Since sin contains only odd powers of Hj we have: 

sm(ojHjT) = Hj sin(wT), (3.4) 
while, since all even powers of Hj equal Hj we can write: 

cas(uHjT) = 1 - Hf[-1 + 1 + \uj 2 T 2 - ...} 

2 

= 1- Hj +H?cos(uT). (3.5) 
We now choose for T a value such that cos(wT) = 0, sin(wT) = 1, getting: 

expi-iuHjT) = 1 - H] - iH 3 . (3.6) 
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The last equation implies that 

exp(-iwHjT) |0j-> <g> |a > = [1 - ff| - iff,-] 1^) ® |a ) 

= |0j)®K-), (3.7) 

as desired. 

We remark now that [iij-,.??/.] = and Hk\4>j) ® |oo) = for A; 5^ J. Accordingly, if consideration is 
given to the hamiltonian i? = 53j=o we nave: 

exp(-iujHT) \(j>j) ® |o ) = |<fo) ® 1%), Vj. (3.8) 

Therefore, we have explicitly exhibited an hamiltonian which performs our game, i.e., it leads to the 
desired purification of our statistical mixture. 

Actually, the purification procedure becomes interesting when one can prepare a desired mixture 
among all those of an equivalence class at-a-distance, as appropriately stressed by Gisin p|. To reach 
this goal a very small change in our formalism is necessary. Let us identify the states \dj) of our equation 
with the internal eigenstates of a system (e.g. the stationary states of an hydrogen atom). One can 
then add to our hamiltonian a term jPcmj where 7 is an appropriately chosen c- number and Pcm is 
the center-of-mass momentum of the system. The evolution induced by the total hamiltonian implies a 
displacement of the system K of an amount governed by the value of 7, so that in the time interval T 
it is brought arbitrarily far from the space region where its interaction with S took place. In brief, the 
auxiliary system is far apart and one can actually use the pure state to prepare the desired statistical 
ensemble of systems S at-a-distance. 

4 A quantum computational example 

It is interesting to notice that, for most cases of interest in quantum computational theory, the outlined 
procedure can be easily implemented by resorting to elementary logical gates. To this purpose, let us 
suppose that S is a qubit, i.e. a two-level system, and let us denote as |0), |1) the computational basis 
of TL S . The controlled-NOT operator acting on qubit S (taken as the control bit) and on the two- 
dimensional manifold spanned by the computational basis states \olq), \a\) of system K (taken as the 
target bit) induces precisely the transformation: 

|0)(8|ao) => |0)®|ao) 

|l)®|ao) => |l)®|oi>, (4.1) 

which is the desired evolution. In this way, we can purify any statistical ensemble belonging to the 
equivalence class of 

£(p,|0);l-p,|l)), 0<p<l, (4.2) 



by starting with an appropriate superposition analogous to the one of Eq. (2.2). 

Let us now consider an arbitrary equivalence class, different from the previous one and containing the 
ensemble (corresponding to the diagonal form of p) : 

£(q,\x + )\l-q,\x-)), < q < 1, (4.3) 

where |x+) and \x-) are a basis obtained from the computational basis |0), |1) by an appropriate "rota- 
tion" of the system: 

\x+) = R s \0), Is-) = Rs\l). (4.4) 
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The circuit that implements the evolution 



\x+) <2> \a ) 
\x-) <8> |a ) 



\x+) ® \a ), 
\x-) ® |ai), 



(4.5) 



leading to the purification of the ensemble, corresponds to a "rotation" Rg on the control bit, followed 
by a controlled-NOT gate and by an inverse "rotation" R$, as shown in the picture. 




Thus, the Hamiltonian that induces the desired evolution can be identified with a rotation in TL S , a 
controlled-NOT operation in H s <g> H K , and finally a counter-rotation in H s . 

In this way, we have identified the appropriate way to purify any statistical ensemble of the two- 
dimensional system S. Useless to say, our procedure can be easily generalized to systems containing 
several qubits and, more in general, to arbitrary quantum systems. 
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